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Abstract —A model of one of the classes of a free-piston Stirling engine, namely a two-piston engine 
with V-shaped positioning of cylinders, is considered. The differential equations describing the 
dynamics of the engine as a self-oscillating system are derived taking sluggishness of the moving engine 
elements (pistons) and redistribution of the actuating medium volumes in the cylinders into account. 
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Internal combustion engines operated on liquid and gaseous fuel (piston and gas turbine ones) are cur¬ 
rently used most widely in different branches of industry. However, during several recent decades scientific 
and practical interest, both in Russia and abroad, is being paid to the machines (engines included) oper¬ 
ated based on another principle, namely external heating using a Stirling cycle [1—4]. Such engines have 
a number of advantages as compared to internal combustion engines, including circulation of the working 
medium in a closed volume, high efficiency, a low level of noise (the audible range included), and a pos¬ 
sibility to use any heat source. 

The Stirling engine is a heat engine operating by a closed thermodynamic cycle where all thermody¬ 
namic processes of compression and expansion take place at different temperature levels and a flow of the 
working medium is controlled by varying its volume. From a design point of view, the Stirling engine con¬ 
sists of expansion and compression cylinders, a regenerator, and a cooler. It should be noted that the 
majority of Russian and foreign papers devoted to the analysis of the Stirling engines that have been pub¬ 
lished by now give only a thermodynamic analysis of their operation cycle. The investigation of displace¬ 
ment of their moving elements (pistons) has been performed, as a rule, based on the kinetostatic analysis, 
i.e., without taking sluggishness of their moving elements into account. 

The authors of [5, 6] took dynamics of moving elements into account and, based on this, the differen¬ 
tial equations of their motion at different sections of a closed cycle were derived for a free-piston two-cyl¬ 
inder engine with coaxial cylinders. The calculations performed based on such model make it possible to 
introduce substantial quantitative and qualitative corrections into the piston motion laws. 

It follows from the thermodynamic calculations, without dynamics of piston displacement taken into 
account, that during the second (isochronous compression) and the fourth (isochronous expansion) cycle 
sections the pistons in both cylinders move with identical velocities preserving a constant distance between 
one another. The dynamic calculation performed in [5, 6] shows that such motion is impossible. This con¬ 
firms the importance of taking the dynamics of moving elements into account when calculating work pro¬ 
cesses in Stirling engines. 

In the present paper, a dynamic model is calculated for a free-piston V-shaped Stirling engine of simple 
operation (Fig. 1): (7, 3) working pistons; (2) a generator; h is the piston thickness, z = x is the length of 
the inner piston surfaces, / is the length of branch pipes, and 
d is the branch pipe diameter. 

As in [5, 6], the following assumptions were made when 
developing the model: hydraulic, mechanical, and heat 
losses are absent; all processes are reversible; regeneration is 
ideal; a steady state mode of engine an operation is consid¬ 
ered. As distinct from a free-piston engine considered in [5, 

6], the dead zone for such engine type turns out to be non¬ 
zero (Fig. 1). 

An ideal operating cycle of the engine consists of four sec¬ 
tions (moves): curve 1—2 correspond to isothermal compres¬ 
sions; curve 2—3, to isochronous compression; curve 3—4 , to 
isothermal expansion; and curve 4—1 , to isochronous expan¬ 
sion. 
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The X symbol marks the distance from the extreme low position of the pis¬ 
ton in the right cylinder (compression cavity) to its extreme top position 
(Fig. 1), where x is the current displacement of this piston. The Zsymbol marks 
the distance from the extreme low position of the piston in the left cylinder 
(expansion cavity) to its top position, where z of the current displacement of the 
left piston. 

The dynamics equations for the four cycle sections are written as follows. 
In the first cycle section (isothermal expansion), the first piston creates a force 
under the action of the external pressure p 0 that exceeds that formed by the 
working medium pressure. The left piston remains immobile. In this case, the 
current volume of the working medium (with the branch pipe volume taken 
into account) is 

F r (x) = y(y-x- + nd 2 l+ n — tan-. 




In the function x, the current value of V r (x) (1) linearly decreases from V rmax (h/ 2) = SX + ncPl + 


d 


h 




tt — tana , corresponding to the extreme low position of the piston to V min [X- - ] = ncPl + n — tan - . In 


this case, the current value of the pressure p r (x) in the compression cavity increases. If we consider the 
working medium to be an ideal gas, the pressure p r (x) and the current value of the volume V r (x) (1) are 
interrelated with the Mendeleev—Clapeyron equation 


p r (x) V r (x) = G n RT mm , 



where G 0 is the total mass of the working medium, R = R\i/M is the universal gas constant, p is the specific 
gas constant, Mis the molar mass of the working medium, and T min is the working medium temperature 
in the right cylinder. 

Based on Eq. (2), let us determine the expression for the current value of the pressure in the compres¬ 
sion cavity 


Pr(x) 


G,R Tmin 

V r (x) 



Thus, two forces act on the piston in the right cylinder: the external force p 0 S from below, and the force 
p r (x)S (3) from above (outwards the right cylinder). Since these forces are oppositely directed (along the 
cylinder axis), the resultant force acting on the piston is F u = p 0 S — p r (x). Thereby, the differential equa¬ 
tion for the piston motion has the form 


mx 



Go^T min 

V r (x) 



h 

It is necessary to integrate Eq. (4) with the entry conditions: x(0) = X -, x(0) = 0. The duration of 

the first cycle section is determined from the condition of the piston halt x(t { ) = 1. Based on this, the pis¬ 
ton coordinate at the halt moment is determined, x 10 = x(t x ). 

At the second cycle section, both pistons move. The left piston starts moving down (from the regener¬ 
ator), while the right piston moves up. The current value of the working medium volume in the expansion 
cavity is determined with the equation 

v { (z) = s{z-z-^) +7r// r + 7tytan|. (5) 
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It is seen from the cavity volume that expansion decreases. In this case, the working medium extruded 
by the right pistol passes via the regenerator to the left cavity where, as the pistons move, its temperature 
increases according to the dependence obtained in [5] 

Uz) = T min + ^( Zl -z), (6) 

where z x <Z — h/ 2 is the coordinate of the extreme top position that was occupied by the left piston before 
the beginning of the second cycle section; C v is the specific heat of the working medium, and a is the coef¬ 
ficient of proportionality between the amount of the heat and the volume of the working medium passing 
via the regenerator in the process of its flow from the compression cavity (the right cylinder) to the expan¬ 
sion cavity (the left cylinder). According to Eq. (6), when the lift cylinder goes down, the temperature of 
the working medium in the expansion cavity increases. 

The values of V { (z) (5) of the pressure p { (z) and 7](z) (6) in the expansion cavity are interrelated by the 
Mendeleev—Clapeyron equation 


Px{z)V x {z) = G l (z)RT l (z), (7) 

where G { (z) = G 0 — G r (x) is the current mass of the working medium in the expansion cavity and G r (x) is 
the current value of the working medium volume in the right cylinder (the compression cavity). 

In this cycle section, G 0 in the Mendeleev—Clapeyron equation for the right cylinder has to be replaced 
by G r (x) and, as the result, it takes on the form 

p r (x)V r (x) = G r (x)RT min . (8) 


From Eqs. (7) and (8), let us determine the pressures in the left and right cavities 


Pl(z) = 


Gi(z)RT l (z) 

Vi(z) 


Pr(x) = 


G r (x)RT mm 


K(x) 



According to [7], the velocity of the working medium flow from the right cavity to the left one via the 
regenerator is related to the pressure difference in the cavities by the following equation 


Pi(z)-p t (x) = Y 


dG x (z) 

1 , 5 


( 10 ) 


where y is the coefficient of the regenerator resistance. 

The motion equations for the left and the right pistons have the form 

mi = \Pi(z)~P 0 \S, mx = [/?o(z)-yMX>]S- (H) 

Equations (10) and (11) form a closed set of nonlinear differential equations with respect to three func¬ 
tions: the current values of the coordinates of the right and left pistons x(t ), z(t ), and the current value of 
the working medium mass in the left piston G x (x , z, t). 

The entry conditions for this set of equations (the moment of the beginning of the second cycle t x cor¬ 
responding to the end of the first cycle) are specified as follows: G l (x l , z x , t x ) = 0 and z(t x ) = z x , x(t x ) = x x . 

The moment of the end of the second cycle section is determined from the condition 
G { (x 2 , z 2 , h) = ^ 1 (^ 2 ) = G 0 ; i.e., all working medium has overflown to the left cylinder. The corresponding 
piston coordinates z 2 = z(t 2 ), x 2 = x(t 2 ) determine the extreme low position of the left piston and the 
extreme top position of the right piston. In this case, the temperature in the expansion cavity reaches its 
maximum determined, according to Eq. (6), with the formula 

T m ax = 7 \(Z 2 ) = T mm + ^( Zl - Z 2 ). (12) 

The temperature in the cavity of the left cylinder at the third cycle section remains constant and is equal 
to the maximal value T max (12). In this case, the right piston remains immobile (x 2 = x(t 2 )), while the left 
one continues going down owing to an increase in the working medium volume in the process of its heating 
at the constant temperature T max . 

The volume of the working medium V { (z) determined with Eq. (5) and the pressure p x (z) are interrelated 
by the Mendeleev—Clapeyron equation p x (z) V { (z) = G 0 RT max . From here, let us determine the pressure in 
the function z m - P\(z) = GQRT max /V { (z)- 
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The motion equation for the left piston is 

ml = [p { (z)-Po\S. (13) 

Equation (13) has to be integrated with the entry conditions z(t 2 ) = Z 2 and z(t 2 ) = z 2 , where z 2 is the 

velocity of z determined from the solution of set of equations (10) and (11) of the previous cycle section 
calculated at the moment of its termination t 2 . 

The end of the third cycle corresponds to a halt of the left piston z(t 3 ) = 0. From this condition, the 
position of the left piston which it occupies in the end of the third cycle section is determined (z 3 = z(t 3 )). 

The fourth section is calculated in the manner analogous to that for the second section. The difference 
is as follows: the left piston is now going up, while the right one goes down, and the working medium flows 
via the regenerator from the left cylinder to the right one. 

The set of differential equations describing the dynamics of this section is written as 

mx = [p r (x)-p 0 \S, ml = [p 0 ~Pi(z)]S, p { (z)-p,(x) = y dG 'j z \ (14) 

at 

It is necessary to integrate the set of differential equations (14) with the entry conditions x(t 3 ) = x 2 , 
x(t 3 ) = 0, zih) = z 3 , z(h) = 0, and G x (x 2 , z 3 ) = G 0 . 

The moment of the fourth section termination or the period of the entire cycle T is equivalently deter¬ 
mined from the conditions x(T) = x u z(T) = z x , and G l (x l , z x , T) = 0. 

Thus, the cycle dynamics in the first and the third section is described by differential equations (4) and 

(13) , respectively; that in the second and fourth sections is described by differential equation sets (10), (11) 
and (14), respectively. Though differential equations (4) and (13) are integrated in quadratures, the cor¬ 
responding analytical equations are rather complex. As for the differential equation sets (10), (11) and 

(14) , they are not integrated in quadratures. 

Thus, it is necessary to use numerical methods for the analysis of the cycle dynamics. 
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